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Chemo-spectrophotometric evolution of spiral galaxies: 
III. Abundance and colour gradients in discs 
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ABSTRACT 

We study the relations between luminosity and chemical abundance profiles of spi- 
ral galaxies, using detailed models for the chemical and spectro-photometric evolution 
of galactic discs. The models are "calibrated" on the Milky Way disc and are success- 
fully extended to other discs with the help of simple "scaling" relations, obtained in 
the framework of semi-analytic models of galaxy formation. We find that our models 
exhibit oxygen abundance gradients that increase in absolute value with decreasing 
disc luminosity (when expressed in dex/kpc) and are independent of disc luminosity 
(when expressed in dex/scalelength), both in agreement with observations. We notice 
an important strong correlation between abundance gradient and disc scalelength. 
These results support the idea of "homologuous evolution" of galactic discs. 
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1 INTRODUCTION 

Almost all large spiral galaxies present sizeable radial abun- 
dance gradients. This is, for instance, the case for the 
Milky Way disc, showing an oxygen abundance gradient of 
dlog(O/H)/dR~-0.07 dex/kpc in both its gaseous and stel- 
lar components. The origin of these gradients is still a mat- 
ter of debate. A radial variation of the star formation rate 
(SFR), or the existence of radial gas flows, or a combina- 
tion of these processes, can lead to abundance gradients in 
discs (e.g. Lacey and Fall 1985, Koeppen 1994, Edmunds and 
Greenhow 1995, Tsujimoto et al. 1995, Firmani et al. 1996, 
Molla et al. 1997, Matteucci and Chiappini 1998). On the 
other hand, the presence of a central bar inducing large scale 
mixing through radial gas flows tends to level out preexist- 
ing abundance gradients (e.g. Dutil and Roy 1999). Despite 
several studies in the 90ies (e.g. Friedli et al. 1998 and ref- 
erences therein) the relative importance of these processes 
has not been clarified yet. 

Several large spectophotometric surveys and literature 
compilations have appeared in recent years (Vila-Costas and 
Edmunds 1992, Zaritsky et al. 1994, Garnett et al. 1997, 
van Zee et al. 1998). These studies allowed to establish sev- 
eral important features of abundances in spirals (Henry and 
Worthey 1999, and references therein). In particular, it is es- 
tablished now that the absolute abundance at a given galac- 
tocentric distance is correlated to galaxy luminosity (and 
rotational velocity). Also, when abundance gradients are 
expressed in dex/kpc there is a tendency of smaller discs 
to exhibit larger gradients; but when they are expressed in 



dex/Rd (where Rd is the disc exponential scalelength), that 
correlation disappears (Garnett et al. 1997). 

In a series of recent papers (Boissier and Prantzos 
1999a, 1999b; hereafter BP99a and BP99b, respectively) we 
presented a detailed model of the chemical and spectropho- 
tometric evolution of spiral galaxies. Adopting some simple 
prescriptions for the star formation and infall rates, we first 
applied the model to the Milky Way disc and showed that 
it reproduces all major observables of our Galaxy (BP99a), 
including radial profiles of gas, stars, SFR and oxygen abun- 
dances. Assuming that the Milky Way is a typical spiral 
galaxy, we extended then the model to the study of other 
spirals (BP99b) with the help of some simple "scaling laws" , 
which have been established in the framework of the Cold 
Dark Matter scenarios for galaxy formation. We showed 
that this simple "scaled" model reproduces quite satisfac- 
torily most of the observed integrated properties of present 
day spirals: disc sizes and central surface brigthness, Tully- 
Fisher relations in various wavelength bands, colour-colour 
and colour-magnitude relations, gas fractions vs. magnitudes 
and colours, abundances vs. local and integrated properties, 
as well as spectra for different galactic rotational velocities. 
Despite the extremely simple nature of our models, we find 
satisfactory agreement with all those observables, provided 
the timescale for star formation in low mass discs is longer 
than for more massive ones. 

In this third paper of the series we explore the implica- 
tions of our models for the radial profiles of spiral galaxies 
and, in particular, the abundance and photometry gradients. 
The plan of the paper is as follows: In Sec. 2 we present 
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the basic ingredients and the underlying assumptions of the 
model. The results of the models concerning the evolution 
of the radial oxygen abundance profiles and the final abun- 
dance and luminosity gradients are presented in Sec. 3. In 
Sec. 4 we compare our results to observations; in particular, 
we show that there is a clear correlation between abundance 
gradient expressed in dex/kpc and disc scalelength, both ob- 
servationnaly and in our models. The implications of these 
results for disc evolution are discussed in Sec. 5. 



2 THE MODEL 

In this section we recall the main ingredients and the under- 
lying asumptions of our chemical and spectrophotometric 
evolution models. 

The galactic disc is simulated as an ensemble of con- 
centric, independently evolving rings, gradually built up 
by infall of primordial composition. The chemical evolu- 
tion of each zone is followed by solving the appropriate set 
of integro-differential equations (Tinsley 1980), without the 
Instantaneous Recycling Approximation. Stellar yields are 
from Woosley and Weaver (1995) for massive stars and Ren- 
zini and Voli (1981) for intermediate mass stars. Fe produc- 
ing SNIa are included, their rate being calculated with the 
prescription of Matteucci and Greggio (1986). The adopted 
stellar Initial Mass Function (IMF) is a multi-slope power- 
law between 0.1 Mq and 100 M Q from the work of Kroupa 
et al. (1993). 

The spectrophotometric evolution is followed in a self- 
consistent way, i.e. with the SFR *(t) and metallicity Z(t) 
of each zone determined by the chemical evolution, and the 
same IMF. The stellar lifetimes, evolutionary tracks and 
spectra are metallicity dependent; the first two are from the 
Geneva library (Schaller et al. 1992, Charbonnel et al. 1996) 
and the latter from Lejeune et al. (1997). Dust absorption 
is included according to the prescriptions of Guiderdoni et 
al. (1998) and assuming a "sandwich" configuration for the 
stars and dust layers (Calzetti et al. 1994). 

The star formation rate (SFR) is locally given by a 
Schmidt-type law, i.e proportional to some power of the gas 
surface density E 9 and varies with galactocentric radius R 



j0 e 



-R/Rd 



(4) 



*(t,i?) = aZ g (t,Ry- 5 V(R) R~ 



(1) 



where V(R) is the circular velocity at radius R. This radial 
dependence of the SFR is suggested by the theory of star 
formation induced by density waves in spiral galaxies (e.g. 
Wyse and Silk 1989). 

The infall rate is assumed to be exponentially decreas- 
ing in time, i.e. 



f(t,R) = A(R) e 



-t/T(R) 



(2) 



with t(Rq—8 kpc) = 7 Gyr in the solar neighborhood, in or- 
der to reproduce the local G-dwarf metallicity distribution. 
The coefficient A(R) is obtained by the requirement that at 
time T=13.5 Gyr the current mass profile of the disc E(i?) 
is obtained, i.e. 



f(t,R) = E(R) 



(3) 



In the case of the Milky Way disc one has for the central 
surface density Eo=1150 Mq pc -2 and for the disc scale- 
length Rd=2.6 kpc, leading to an overall disc mass Md ~5 
10 10 Mq (see BP99a and references therein). 

The efficiency a of the SFR (Eq. 1) is fixed by the 
requirement that the observed local gas fraction a g (Ro=8 
kpc)~0.2, is reproduced at T=13.5 Gyr. We consider then 
that the really "free" parameters of the model are the radial 
dependence of the infall timescale r(R) and of the SFR. It 
turns out that the number of observables explained by the 
model is much larger than the number of free parameters. In 
particular the model reproduces present day "global" prop- 
erties (profiles of gas, O/H, SFR, and supernova rates), as 
well as the current disc luminosities in various wavelength 
bands and the corresponding radial profiles of gas, stars, 
SFR and metal abundances; moreover, the adopted inside- 
out star forming scheme leads to a scalelength of ~4 kpc in 
the B-band and ~2.6 kpc in the K-band, in agreement with 
observations (see BP99a). 

In order to extend the model to other disc galaxies 
we adopt the "scaling properties" derived by Mo, Mao and 
White (1998, hereafter MMW98) in the framework of the 
Cold Dark Matter (CDM) scenario for galaxy formation. 
According to this scenario, primordial density fluctuations 
give rise to haloes of non-baryonic dark matter of mass M, 
within which baryonic gas condenses later and forms discs of 
maximum circular velocity Vc- It turns out that disc profiles 
can be expressed in terms of only two parameters: rotational 
velocity Vc (measuring the mass of the halo and, by assum- 
ing a constant halo/disc mass ratio, also the mass of the 
disc) and spin parameter A (measuring the specific angular 
momentum of the halo). In fact, a third parameter, the red- 
shift of galaxy formation (depending on galaxy's mass) is 
playing a key role in all hierarchical models of galaxy for- 
mation. However, since the term "time of galaxy formation" 
is not well defined (is it the time that the first generation 
of stars form? or the time that some fraction of the stars, 
e.g. 50% form?) we prefer to ignore it and assume that all 
discs start forming their stars at the same time, but not at 
the same rate. In that case the profile of a given disc can be 
expressed in terms of the one of our Galaxy (the parameters 
of which are designated hereafter by index G): 
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with 



where we have adopted Ag=0.03 (see BP99b). Numerical 
simulations show that the A distribution peaks around the 
value 0.04 (e.g. MMW98). The absolute value of A G is of 
little importance as far as it is close to that peak, since our 
results depend only on the ratio A/Ag- 

Eqs. 5 and 6 allow to describe the mass profile of a 
galactic disc in terms of the one of our Galaxy and of two 
parameters: Vc and A. The range of observed values for the 
former parameter is 80-360 km/s, whereas for the latter nu- 
merical simulations give values in the O.OTO.l range, the 
distribution peaking around A ~ 0.04 (MMW98). Although 
it is not clear yet whether Vc and A are independent quan- 
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Figure 1. Evolution of oxygen abundance profiles in our models, at times t=3 Gyr (dotted curves), 7.5 Gyr (dash-dotted curves) and 
13.5 Gyr (solid curves). Models are arranged according to values of maximum rotational velocity Vc (increasing from top to bottom) 
and spin parameter A (increasing from left to right). In the case of the Milky Way model (Vc=220 km/s, A=0.03) observations of 
oxygen abundances in HII regions, planetary nebulae and B-stars are also plotted; the observed metallicity gradient dlog(0/H)/dR~- 
0.07 dex/kpc is well reproduced. Notice that the absolute values of A do not enter directly the calculations, since the model discs are 
scaled to the one of the Milky Way through the relations (5) and (6); in other terms, if the Milky Way A is taken to be 0.06, the same 
results would be obtained for discs with A = 0.02, 0.06, 0.10, 0.14 and 0.18, respectively. 



tities, we treat them here as such and construct a grid of 
25 models caracterised by Vc = 80, 150, 220, 290, 360 km/s 
and A = 0.01, 0.03, 0.05, 0.07, 0.09, respectively. [Notice: if 
Ag=0.06 is adopted for the Milky Way, our model results 
would be the same, but they would correspond to values of 
A twice as large, i.e. 0.02, 0.06, 0.10, 0.14 and 0.18, respec- 
tively]. Increasing values of Vc correspond to more massive 
discs and increasing values of A to more extended discs. 

As discussed in BP99b, the resulting disc radii and cen- 
tral surface brightness are in excellent agreement with ob- 
servations, except for the case of A=0.01. This "unphysical" 
value leads to galaxies ressembling to bulges or ellipticals, 
rather than discs. We prefer, however, to keep this A value 
in our grid, for illustration purposes. 

The two main ingredients of the model, namely the Star 
Formation Rate "if(R) and the infall time-scale t(R), are af- 
fected by the adopted scaling of disc properties in the fol- 
lowing way: 

For the SFR we adopt the prescription of Eq. (1), with 
the same efficiency a as in the case of the Milky Way (i.e. 
the SFR is not a free parameter of the model). In order 
to have an accurate evaluation of V(R) across the disc, we 
calculate it as the sum of the contributions of the disc (with 
the surface density profile of Eq. 4) and of the dark halo, 



with a volume density profile of a non-singular isothermal 
sphere (see BP99b). 

The infall time scale is assumed to increase with both 
surface density (i.e. the denser inner zones are formed more 
rapidly) and with galaxy's mass, i.e. r[Md, £(/?)]. In both 
cases it is the larger gravitational potential that induces a 
more rapid infall. The radial dependence of r on S(i?) is 
calibrated on the Milky Way, while the mass dependence of 
r is adjusted as to reproduce several of the properties of the 
galactic discs (see BP99b). The adopted prescription allows 
to keep some simple scaling relations for the infall in our 
models and the number of free parameters as small as possi- 
ble. The fact that the adopted simple prescription provides 
a satisfactory agreement with several observed relationships 
in spirals (see BP99b) makes us feel that it could be ulti- 
mately justified by theory or numerical simulations of disc 
formation. 

We stress at this point that in our models we ignore 
the possibility of radial inflows in gaseous discs, resulting 
e.g. by viscosity or by infalling gas with specific angular mo- 
mentum different from the one of the underlying disc; in 
both cases, the resulting redistribution of angular momen- 
tum leads to radial mass flows. The magnitude of the effect 
is difficult to evaluate, because of our poor understanding 
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of viscosity and our ignorance of the kinematics of the in- 
falling gas. Models with radial inflows have been explored 
in the past (Mayor and Vigroux 1981; Lacey and Fall 1986; 
Clarke 1989; Chamcham and Tayler 1994); at the present 
stage of our knowledge, introduction of radial inflows in the 
models would imply even more free parameters and make 
impossible the study of a radial variation in the efficiency of 
the SFR. Since it is well known that bars induce radial mix- 
ing and reduce the abundance gradients (e.g. Dutil and Roy 
1999), we do not consider bared galaxies in what follows. 



3 RESULTS AND COMPARISON TO 
OBSERVATIONS 

3.1 Model results 

The evolution of the oxygen abundance profiles in our mod- 
els is presented in Fig. 1. Profiles are shown at three ages: 
T=3, 7.5 and 13.5 Gyr, respectively, covering most of the 
galaxian history. It should be noticed that the observed oxy- 
gen profile in the Milky Way disc is nicely reproduced by the 
corresponding model (Vc=220 km/s, A=0.03), as analysed 
in detail in BP99a. The most important features of Fig. 1 
are the following: 

1) Because of the inside-out star forming scheme 
adopted in our work, the abundance profile is quite steep 
early on (when stars and metals are formed only in the inner 
disc) and flattens with time (as star formation "migrates" 
progressively outwards) . 

2) The final abundance gradients (expressed in 
dex/kpc) are, in general, larger for the lower mass discs. The 
reason is that, in such small galaxies, the S s /R factor in 
the adopted SFR can create large differences in abundances 
between regions relatively close to each other. Indeed, be- 
cause of the small scalelength of these galaxies, the gas sur- 
face density varies a lot within a short distance, and so does 
the corresponding SFR. On the contrary, for massive and 
large discs, the gas surface density and corresponding SFR 
vary slowly within one kpc and produce flatter abundance 
gradients (in dex/kpc). 

3) The more massive the galaxy, the higher is the cen- 
tral (and average) abundance. This is due to the V/R factor 
of the adopted SFR, making massive discs more efficient 
in forming stars and metals at a given galactocentric dis- 
tance. This factor explains also why the most compact discs 
(A=0.01) have very high abundances. However, as analysed 
in BP99b, that case does not produce "physical" results, 
since the obtained galaxies ressemble more to ellipticals or 
bulges than to spirals. We show the results for A=0.01 in 
Fig. 1 for illustration purposes, but in the following we shall 
discuss only results for the other cases. In BP99b we show 
that this increase of metallicity with disc mass and luminos- 
ity compares well to the data of Zaritsky et al. (1994). 

4) The most efficient galaxies in forming stars and met- 
als are the compact (low A) and massive (large Vc) ones. In 
their innermost regions, metallicity tends to saturate, to a 
value of log(0/H) + 12~9.4. This is a well known feature of 
galactic chemical evolution models taking into account the 
finite lifetime of stars (see Prantzos and Aubert 1995). At 
late times, when most of the gas is exhausted and little star 
formation or metal production takes place, large amounts 
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Figure 2. Results of our models at an age T=13.5 Gyr. They 
are plotted as a function of maximal rotational velocity Vq and 
spin parameter A (0.03 solid curves, 0.05 dotted curves, 0.07 short- 
dashed curves and 0.09 long-dashed curves). From top to bottom: 
oxygen abundance gradient (in dex/kpc); absolute B-magnitude 
Mb; disc scalelength in the B-band Rg; radius where B-band 
surface brightness is 25 mag arcsec -2 (R25); colour gradient B- 
K (in mag/kpc). In all panels, the symbol indicates values 
corresponding to the Milky Way model. 



of metal-poor gas are returned to the interstellar medium; 
they are ejected by the numerous low-mass and long-lived 
stars that were abundantly formed in the early days of the 
galaxy. This metal-poor gas dilutes the metal abundances 
and leads to a "saturation" of their value. This feature can- 
not be revealed by models using the Instantaneous Recycling 
Approximation. For that reason, the oxygen profiles of the 
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Figure 3. Disc scalclcngths in various wavelength bands (left) and disc scalelength ratios (middle and right) as a function of B-magnitudc. 
Our results (curves, parametrised with values of the spin parameter A as indicated on the left) are compared to the observations of de Jong 
(1996). Model scalelengths in the middle panels are calculated with extinction taken into account, while in the right panels extinction is 
neglected. In all panels, the filled circle indicates results for the Milky Way disc. 



A=0.01 case are fiat, since star formation is very rapid and 
efficient all over the (small) disc. 

Figure 2 summarizes a number of important features of 
our models at an age of T=13.5 Gyr. 

a) The oxygen abundance gradient (in dex/kpc) is a 
monotonically increasing function of Vc- Notice that gra- 
dients are negative and more massive discs have smaller 
absolute values of d(log(0/H))/dR, i.e. flatter abundance 
profiles. Also, abundance gradients depend on the A value, 
especially in low mass discs: more compact discs have more 
important (in absolute value) gradients, while extended discs 
have flat abundance profiles, fn both cases, the fact that 
smaller discs have larger abundance gradients (in dex/kpc) 
is due to the E /R factor in the adopted star formation 
rate, as explained in point 2 above. 

b) The B-magnitude depends mainly on the mass of the 
galaxy (i.e. rotational velocity Vc) and very little on the spin 
parameter A. As analysed in detail in BP99b, a Tully-Fisher 
relation with a slope similar to the observed one in field 
spirals and small dispersion is obtained in our models. This 
is a success of semi-analytical models, since the TF relation 
is "built-in" as a boundary condition (the disc mass being 
proportional to Vc)- 

c) The final B-band scalelength is a monotonically in- 
creasing function of Vc and A. This again is due to the 
"boundary conditions" imposed by the semi-analytic models 
adopted in our work (see Sec. 2). 

d) R25, the radius where B-band surface brightness is 
25 mag/arcsec 2 , is also found to increase monotonically with 
Vc and A. 



e) Finally, discs with A > 0.05 have small colour gra- 
dients, independently of their mass (lower panel in Fig. 2). 
This is because in those galaxies the whole disc is formed 
relatively late (in timescales >5 Gyr for Vc <220 km/s, see 
Fig. 5 in BP99b) and stellar populations in the inner disc 
are not much older than in the outer regions. 



3.2 Comparison to observations of luminosity 
profiles 

The exponential scalelength of the Milky Way disc was a 
subject of controversy for several years. Indeed, observa- 
tions in different wavelength bands lead to different results: a 
value of Rs ~4-5 kpc is found in the B-band, while smaller 
values are found in longer wavelengths (Rj ~2.5 kpc). In 
BP99a we showed that these different values result natu- 
rally from schemes of "inside-out" formation of the disc, 
that manage to reproduce the observed oxygen abundance 
gradient of dlog(O/H)/dR~-0.07 dex/kpc. In fact, light in 
the I- or K- bands traces closely the old stellar population, 
which dominates the mass of the disc in all radii. Light in 
the B-band traces a younger stellar population, formed in 
the past 1-2 Gyr from a gas with a surface density profile 
flatter than the one of the old population. Thus, colour and 
abundance gradients are found to be naturally correlated in 
the case of the Milky Way disc. 

In this section we explore the results of our models con- 
cerning the disc scalelengths in various wavelength bands. 
Our aim is: 1) to see whether our results are compatible 
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with the observations and 2) to check whether there are any 
similarities with the case of the Milky Way. 

On the left part of Fig. 3 we present the disc scalelength 
of our models in the B-, V- and K- bands (evaluated with the 
effects of extinction taken into account) and compare them 
to the observations of De Jong (1996). De Jong has decom- 
posed the luminosity profiles of his ~80 face-on galaxies in 
bulge and disc components and we plot in Fig. 3 only the 
corresponding disc component, directly comparable to our 
results. It can be seen that the model disc scalelengths are in 
excellent agreement with the observed ones as a function of 
Ms . More massive and bright discs have larger scalelengths 
on the average, while the variation of the spin parameter A 
allows to reproduce the observed scatter (perhaps not quite 
efficiently in the B-band). 

On the right part of Fig. 3 we plot the disc scalelength 
ratios as a function of Ms. Since extinction may affect the 
results, we plot our results with extinction included (mid- 
dle panels) and neglected (right panels). It can be seen that 
the observed Rs/Rr- and Rv/Rk rates are on the average 
larger than 1 (i.e. there is a colour gradient); also, there is 
a marginal trend of increasing ratios with galaxy luminos- 
ity. Our results are in agreement with the observed trends, 
in particular when extinction is included. We find that the 
case of A=0.03 leads to the most extreme scalelength ratios; 
other A values (i.e. less compact discs) produce more mod- 
est colour gradients. In that respect then, our Milky Way 
is a rather compact disc having larger than average colour 
gradients. 



3.3 Comparison to observations of abundance 
gradients 

In this section we compare our results to the observations of 
Garnett et al. (1997) and van Zee et al. (1998) concerning 
abundance gradients in external spirals (i.e. observations of 
O/H in HII regions as a function of galactocentric distance). 

Garnett et al. (1997) combine data obtained in several 
previous works (Vila-Costa and Edmunds 1992, Zaritsky et 
al. 1994 and Ryder 1995) and make several adjustements 
in order to introduce a greater degree of uniformity. From 
this final sample we selected a sub-sample of 17 non-barred 
galaxies, since those compare better to our models (bars can 
induce mixing of gas along the disc and flattening of abun- 
dance gradients in relatively short timescales, e.g. Friedli 
et al. 1998 and references therein; our models do not take 
such effects into account). On the other hand, van Zee et 
al. (1998) combine data from the literature (including their 
own observations for low metallicity HII regions) and derive 
abundance gradients for 11 more spirals. 

Our results are compared to observations in Fig. 4, 
where abundance gradients are plotted as a function of in- 
tegrated B-magnitude. 

In the upper panel, abundance gradients are expressed 
in dex/kpc. Observations show that luminous discs have 
small gradients, while as we go to low luminosity ones the 
absolute value of the gradients and their dispersion increase. 
This trend is fairly well reproduced by our models. The rea- 
son can be understood by comparing to the upper panel of 
Fig. 2, where our results are plotted as a function of Vc 
(which correlates well to Ms through the Tully-Fisher rela- 
tion, BP99b). As analysed in Sec. 3.1, the Eg' 5 /R factor in 
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Figure 4. Oxygen abundance gradients of our models at T=13.5 
Gyr (open symbols, corresponding to the values of the spin pa- 
rameter A as indicated in the upper panel and connected with 
solid curves). They are plotted as a function of the B-magnitudc 
and compared to observations (filled symbols, with appropriate 
references given in each panel). Upper panel: Oxygen gradient in 
dex/kpc; Middle panel: Oxygen gradient in dex/Rs ; Lower panel: 
Oxygen gradient in dex/R,25. 



the adopted SFR creates large abundance variations within 
short distances in small galaxies; no important gradients can 
be created in massive discs, where neighboring regions differ 
little in E 9 . As for the increased dispersion at low luminosi- 
ties (and low Vc), it is due to the effect of the A parameter: 
for large A values, even small mass discs are quite extended 
(Eq. 5 in Sec. 2) and the E 1,5 /R factor becomes less and 
less efficient in creating abundance gradients (since £ varies 
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little as a function of R). Thus, in low mass discs, both the 
mass and the spin parameter A drive the abundance gradient. 

As stressed in BP99b, discs with A > 0.08 correspond 
to Low Surface Brightness galaxies (LSBs, see also Jimenez 
et al. 1998). According to DeBlok and van der Hulst (1998) 
LSBs have negligible abundance gradients. Our results of 
Fig. 4 are quite encouraging in that respect: indeed, had we 
run models with A > 0.09 we would expect to find smaller 
and smaller abundance gradients at all galaxy luminosities. 
We postpone, however, the detailed study of LSBs in a forth- 
coming paper. 

In the middle panel of Fig. 4, the oxygen gradients are 
expressed in dex/Rs. When expressed in this unit, the ob- 
served abundance gradients show no more any correlation to 
Mb, as already noticed in Garnett et al. (1997). Moreover, a 
considerable dispersion is obtained for all Ms values, while 
the average gradient is ~-0.2 dex/Rs. Our results (values 
in the upper panel of Fig. 2, multiplied by the correspond- 
ing Rb values of the middle panel of Fig. 2) show also no 
correlation with Ms, in fair agreement with observations. 
However, the average value is slightly larger than observed 
(~-0.25 dex/Rs). Since the estimates of Rs, both in obser- 
vations and in our models, may be affected by extinction, 
we consider that the agreement with the data is quite sat- 
isfactory. Another difficulty may stem from the fact that 
relatively few discs can be fit with perfect exponentials; ac- 
cording to Courteau and Rix (1999) this happens for only 
~20% of the disc galaxies. Finally, we obtain no significant 
dispersion in our models, contrary to obervations. We sus- 
pect that the difficulties to evaluate the effects of extinction 
or to fit the profiles with a perfect exponential may be (at 
least partly) responsible for the observed scatter. 

In the lower panel of Fig. 4 abundance gradients are ex- 
pressed in dex/R25. Again, observations show no trend with 
Ms and a considerable dispersion. Our results also show no 
correlation with Ms, while the average model value (~-0.8 
dex/R25) is in perfect agreement with observations. This is 
most encouraging, since R25 is less affected by considerations 
on dust extinction or fit to an exponential profile. Indeed, we 
think that the obtained average value of dlog(0/H)/dR=- 
0.8 dex/R25 is a more robust result than the -0.25 dex/Rs 
value obtained in the previous paragraph. However, the ob- 
served dispersion is not reproduced by our models; only part 
of it can be attributed to the effects of the spin parameter 
A. 



3.4 Abundance and colour gradients vs. 
scalelength 

In the previous section it was shown that: 

a) abundance gradients (in dex/kpc) become more im- 
portant and present a larger dispersion in low luminosity 
discs. 

b) Abundance gradients in dex/Rs are independent of 
disc luminosity, and so is the corresponding dispersion. 

Points (a) and (b) immediately suggest that there must 
be a one-to-one correlation between abundance gradients 
and scalelength Rs- The correlation must be relatively tight, 
since the observed dispersion (middle panel of Fig. 4) is rela- 
tively small. Although this conclusion is an immediate con- 
sequence of the observations in Fig. 4, it has never been 
stated explicitly (to our knowledge). For the same reason, 
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Figure 5. Predictions of our models for the gradients of metallic- 
ity (in dex/kpc) vs 1/Rs (left panels) and I/R25 (right panels). 
In the upper panels our models are parametrised by Vc values 
and in the lower panels by A values. Data are from Garnett et al. 
(1997), van Zee et al. (1998) and Zaritsky et al. (1994). A clear 
correlation is found between abundance gradient (in dex/kpc) and 
1/Rs- 



a one-to-one correlation must exist between abundance gra- 
dients and R25. However, the observed scatter (lower panel 
in Fig. 4) shows that one should not expect a very tight 
correlation in that case. 

In Fig. 5 we plot the oxygen abundance gradients (in 
dex/kpc) as a function of 1/Rs (which is proportional to 
a "magnitude gradient" in mag/kpc) and of I/R25, both 
from observations and from our models. In the upper pan- 
els our results are parametrised by the values of rotational 
velocity Vc, while in the lower panels they are parametrised 
by the values of the spin parameter A. It can be seen that 
our models show an excellent correlation between the abun- 
dance gradients and 1/Rs: smaller discs have larger gradi- 
ents and the results depend little on A. In principle, knowing 
the B-band scalelength of a disc, one should be able to deter- 
mine the corresponding abundance gradient and vice- versa! 
In practice, however, observations show a scatter around the 
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Figure 6. Predictions of our models at T=13.5 Gyr for the re- 
lation between metallicity gradient (in dex/kpc) and colour gra- 
dient B-K (in mag/kpc). The grid of our models is parametrised 
by V c (360, 290, 220, 150 and 80 km/s, from top to bottom) 
and A (0.03, 0.05, 0.07 and 0.09, from left to right). Massive and 
extended discs have negligible metallicity and colour gradients. 



observed trend. Moreover, our results seem to agree better 
with the data of Garnett et al. (1997) than with those of 
van Zee et al. (1998). The latter show a less steep depen- 
dence of the gradient on 1 /Kb ■ Since large discs "cluster" 
on the upper left corner of the diagram, we suggest that ob- 
servations of the abundance gradients in small discs would 
help to establish the exact form of the correlation. Finally, 
we notice that in our models we find similar correlations 
between abundance gradients and disc scalelengths in other 
wavelength bands (as can be expected from the fact that 
scalelength ratios are close to unity, Fig. 3). 

The situation is not as promising when abundance gra- 
dients are plotted as a function of I/R25. Our models show 
again a correlation, but this time the different A values lead 
to increasingly different abundance gradients as one goes to 
smaller discs. Although the data are compatible with the 
obtained relation, the large observational uncertainties do 
not allow to conclude. Again, further observations for small 
discs are required to clarify the situation. 

Finally, our models predict that, in general, colour gra- 
dients are developped in discs, because of the inside-out star 
formation scheme. For instance, in the case of the Milky Way 
we obtain Kb ~4 kpc and Ra _ ~2.6 kpc, in agreement with 
observations (see BP99a and references therein). As shown 
in Fig. 2 (bottom panel) and discussed in Sec. 3.1, colour 
gradients are rather small for discs with A >0.05; only for 
low A and Vc values do they become important. In Fig. 
6 we present the obtained relationship between abundance 
and B-K gradient, expressed in dex/kpc and mag/kpc, re- 
spectively. Although the resulting grid covers a rather large 
part of the figure, it nevertheless excludes the co-existence 
of large colour and small abundance gradients. Combined 
chemical and photometric analysis of disc profiles should al- 



low to find out whether the proposed scheme is realistic or 
not. 



4 SUMMARY 

In this work, we study the photometric and chemical abun- 
dance profiles resulting from our models of the evolution of 
galactic discs. Our galaxy evolution models are calibrated on 
the Milky Way (BP99a) and reproduce a large body of ob- 
servational data concerning integrated properties of external 
spirals (BP99b) . Here it is shown that they also produce re- 
sults in excellent agreement with observations of abundance 
gradients in spiral galaxies. More specifically: 

1) Because of the inside-out star formation scheme, our 
models predict that abundance and photometric profiles flat- 
ten with time. 

2) Model abundance gradients (in dex/kpc) are more 
important and present a larger dispersion for low luminos- 
ity discs, in fair agreement with observations. We claim that 
the spin parameter A (affecting the size, but not the mass of 
the disc in our models) is responsible for the scatter. Discs 
with larger A values than those studied here should appear 
as Low Surface Brightness galaxies and should exhibit even 
smaller abundance gradients than those displayed in Fig. 3. 
Observations show indeed that LSBs have negligible abun- 
dance gradients (DeBlok and van der Hulst, 1998). 

3) When expressed in dex/Rs and dex/R25, observed 
abundance gradients show no dependence on galaxy lumi- 
nosity; they do show an important scatter, especially in the 
latter case. Our models also show no trend with luminosity, 
while the average values are in fair agreement with obser- 
vations. We fail, however, to reproduce the observed scatter 
(variations induced by the different A values are not suffi- 
cient for that). 

4) The observed and theoretical relation between galaxy 
luminosity and abundance gradients, suggest that there 
must also exist a correlation between the latter and 1/Rs 
or I/R25. Our models show indeed such a correlation: abun- 
dance gradients (in dex/kpc) correlate much better to 1/Rs 
than to Mb. Observations seem to coroborate this idea, al- 
though more data are needed to establish the precise form 
of that correlation. 

5) Our models show no correlation between abundance 
and colour gradients (expressed in dex/kpc and mag/kpc, 
respectively), except for the most compact discs (those with 
A=0.03). Compact discs of low luminosity are found to have 
large abundance and colour gradients. 

The most interesting feature of this work is the anti- 
correlation found between abundance gradients and scale- 
length. Discs with small scalelengths have large abundance 
gradients, independently of their mass. It is the first time (to 
our knowledge) that this relationship is noticed. It should be 
interesting to establish it through further observations, con- 
cerning in particular discs of small scalelengths. Our models 
seem to reproduce that relationship quite naturally (Fig. 5). 

We wish to stress here that our simple chemical evolu- 
tion models treat galactic discs as ensembles of concentric, 
independently evolving rings, i.e. no radial flows ar consid- 
ered. If radial inflows are indeed playing a key role in shaping 
abundance gradients (as sometimes claimed), how is then 
the success of our models to be understood? 
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As explained in Sec. 2, our models are "calibrated" on 
the Milky way disc. In particular, the radial dependence of 
the adopted infall rate is such that, when combined to the 
adopted SFR, reproduces the current profiles of gas, stars 
and metal abundances. The fact that we keep the same pre- 
scriptions for the SFR and infall rate in the other disc mod- 
els, leads obviously to a kind of "homologuous evolution" , as 
already suggested in Garnett et al. (1997). In other words, 
prescriptions adopted for the Milky Way are found to be 
equally successful when applied to other discs. It should be 
interesting to see whether hydrodynamical simulations (like 
e.g. those of Samland et al. 1997, or Contardo et al. 1998) 
manage to produce this kind of "homologuous evolution" in 
galactic discs. We notice, however, that the idea of a "homo- 
loguous evolution" merely states a fact, but does not itself 
explain the origin of the observed gradients. 
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